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INVERSION FORMULAS FOR THE SPHERICAL 
MEANS IN CONSTANT CURVATURE SPACES 

YURI A. ANTIPOV, RICARDO ESTRADA, AND BORIS RUBIN 



Abstract. The work develops further the theory of the following 
^^' inversion problem, which plays the central role in the rapidly de- 

.^fH , veloping area of thermoacoustic tomography and has intimate con- 

nections with PDEs and integral geometry: Reconstruct a function 
^^ . / supported in an n-dimensional ball B , if the spherical means of 

f are known over all geodesic spheres centered on the boundary of 
.^ . B. We propose a new unified approach based on the idea of ana- 

r \ ' lytic continuation. This approach gives explicit inversion formulas 

• . not only for the Euclidean space M" (as in the original set-up) 

"^ I but also for arbitrary constant curvature space X, including the 

C^ , n-dimensional sphere and the hyperbolic space. The results are ap- 

plied to inverse problems for a large class of Euler-Poisson-Darboux 
equations in constant curvature spaces of arbitrary dimension. 

(N' 

> 

(N. 

0\ , 1. Introduction 

a^ 

•/^ \ The paper deals with the spherical mean operator, which is also 

t^^ ■ known as the spherical mean Radon transform. Importance of this 

^ . transformation in analysis and geometry and many of its properties 

(which are still surprising!) were indicated by many authors; see, e.g., 
[11, p. 699], [31, 57]. In recent years an interest to this object has 
grown tremendously in view of a series of challenging problems. One 
of them is characterization of sets of injectivity (and non-injectivity) of 
5^ \ this transform; see [15, 4, 5, 56] and references therein. Another source 

of mathematical problems related to the spherical means is the rapidly 
developing thermoacoustic tomography (TAT), the revolutionary role 
of which in medical imaging was pointed out in many publications; see 
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[l]-[3], [18, 19], [20]-[23], [30, 35], [34]-[36], [38], [46]-[50], [59, 61]. The 
present investigation belongs to this area. 

Setting of the problem and motivation. Let / be an infinitely 
differentiable function with compact support in the open ball B = {x & 
R"^ : |a;| < R}; OB is the boundary of B. We consider the spherical 
mean Radon transform Mf which integrates / over spheres centered 
on dB: 

(1.1) (Mm,t) = ^f i\i-te)de, 

where 

^edB, tGM+ = (0, oo), 

S*"^"*^ is the unit sphere in M" with the area o"„_i, and da stands for the 
usual Lebesgue measure on S"'~^. For the classical Radon transforms, 
their modifications, and applications see, e.g., [12, 14, 26, 13, 29, 42]. 

The general problem of reconstructing / from known data {Mf){^, t) 
on the cylinder dB x K+ is an immediate consequence of the following 
commonly accepted mathematical model of TAT in K^ (see, e.g., [36, 
59, 30] and references therein): 

Given a function c{x), the speed of the ultrasound propagation in the 
tissue, and a function g{^, t), the measured value of the pressure at the 
time t at the transducers location C, & S"^, find a function f{x), the 
initial pressure distribution p{x, 0) (the TAT image), if 

{Ptt = c^{x) Ap for all t > 0, x e M^, 

p{x, 0) = f{x), pt{x, 0) = for all x G R^ 
P{U) = giU) for all ^eS^c m3), t > 0. 

Here, pt and pu are the first and second time derivatives, and A is the 
Laplace operator with respect to the spatial variable x. 

This problem admits immediate generalization to arbitrary dimen- 
sions, more general Riemannian spaces, and a broad class of differential 
equations of the Euler-Poisson-Darboux (EPD) type. In the Euclidean 
case, a thorough discussion of main inversion methods in terms of their 
assumptions and computational features can be found in [36, 59]. In 
the important particular case of constant speed c(x), solution to this 
problem is equivalent to reconstruction of / from its spherical mean 
(1.1). 

Explicit inversion formulas for Mf are of particular interest. For n 
odd, such formulas were obtained by Finch, Patch, and Rakesh in [19] . 
Another derivation was suggested by Palamodov [47, Section 7.5]; see 
also [21, 22, 61]. The corresponding formulas for n even were obtained 
by Finch, Haltmeier, and Rakesh in [18]. An explicit inversion formula 
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which rehes on completely different ideas and covers both odd and even 
cases, was suggested by Kunyansky [38]. 

In spite of the elegance and ingenuity, the derivation of the existing 
inversion formulas for Mf is pretty involved, and basic ideas behind 
it remain mysterious. In view of practical importance, it would be 
desirable to find an independent simple proof of known formulas and 
thus check their correctness. Moreover, the prospective new method 
should be applicable to more general geometric and analytic settings 
and thus lead to further progress. 

A simple proof for n odd was suggested by the third author [54], 
who suggested to treat M as a member of a certain analytic family of 
operators and applied the results to the inverse problem of type (1.2) 
for the more general EPD equation in the case c{x) = const. 

In the present article we suggest a new approach, which is concep- 
tually simple and leads to inversion formulas for Mf for all n > 2. As 
in [54], the key idea is analytic continuation, however, the reasoning is 
different. We extend our method to the similar problem for spherical 
means on the n-dimensional sphere and the hyperbolic space, where 
the theory of EPD equations is also well-developed. This extension 
seems to be new and paves the way to diverse settings, when the rele- 
vant geodesic balls and spherical means are considered in more general 
Riemannian spaces. 

Regarding generalizations to general Riemannian spaces, some com- 
ments are in order. The corresponding wave equations and their EPD 
generalizations were studied in [39, 32, 33], [43]-[45]. For example, the 
wave equation on the n-dimensional sphere S"' has the form [39] 

(1.3) S^u = u^^ + l^^-^j u, (x,w) G 5" X (0,7r), 

where S^ denotes the Beltrami-Laplace operator. The Cauchy problem 
for the relevant EPD equation 

(1.4) D„u = 0, u{x,0) = f{x), u^{x,0) = 0, 
where 

(1.5) DaU = 5xU — Uu)io — (n — 1 + 2a) cotuu^^ + a{n — 1 + a)u^ 

and various modifications were discussed in [9, 10, 24, 32, 33, 44]. The 
problem (1.4) for the case a = 0, corresponding to the usual Darboux 
equation, was studied by Olevskii [43] and also by Kipriyanov and 
Ivanov [32]. Our definition of the EPD-equation on S*" differs from 
that in [32] and agrees with [44]. The particular case a = (1 — n)/2. 
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corresponding to the wave equation (1.3), can be regarded as the spher- 
ical analogue of the TAT model (1.2) with constant speed. 

Plan of the paper and main results. Section 2 contains prelim- 
inaries. Here the main statement is Lemma 2.2. For convenience of 
the reader and better treatment of the subject, we supply this lemma 
with alternative proofs, which are based on different ideas and use dif- 
ferent tools, while leading to the same result. All these are presented 
in Appendix. Section 3 contains derivation of inversion formulas for 
Mf in the Euclidean case. The main inversion results are presented 
in Theorems 3.4 and 3.7; see also modified inversion formulas (3.21), 
(3.22). The results of Section 3 are applied in Section 4 to the Cauchy 
problem for the Euler-Poisson-Darboux equation 

n-\-2a — l 
(1.6) DaU = Au-utt ut = 0, u{x,0) = f{x), Ut{x, 0)=0, 

where / is a smooth function with compact support in the ball B. Using 
the results of Section 3 combined with known properties of Erdelyi- 
Kober fractional integrals, we give explicit solution (Theorem 4.1) to 
the following inverse problem: 

Given the trace u (^,t) of the solution of (1-6) for all {C,,t) on the 
cylindrical surface dB x ]R_|_, reconstruct f{x). 

The particular case a = (1 — n)/2 gives explicit solution to the TAT 
problem (1.2) with constant speed c{x) = 1. 

The spherical mean Radon transform on the n-dimensional unit 
sphere S*" in M""*"^ is studied in Section 5. Inversion formulas for this 
transform are given in Theorems 5.3 and 5.5. The relevant inverse 
problem for the EPD equation on S*" is solved in Section 6. Section 7 
contains derivation of inversion formulas for the spherical mean Radon 
transform in the n-dimensional hyperbolic space. Here the main results 
are given by Theorems 7.3 and 7.5 . 

Acknowledgements. The third author is grateful to Mark Agra- 
novsky, who encouraged him to study this problem, and also to Pe- 
ter Kuchment and Leonid Kunyansky for useful discussions. Special 
thanks go to David Finch, who shared with us his knowledge of the 
subject. 

2. Auxiliary statements 

Notation. We use abbreviation a.c. to denote analytic continuation; 
(T„_i = 27r"/^/r(n/2) is the area of the unit sphere S*""-*^ in M". We 
write d9 (d^) for the usual Lebesgue measure on S*""^ (on dB, resp.); 
[a] denotes the integer part of a real number a; (■)^ means {■)^ if the 
expression in parentheses is positive and zero, otherwise. 
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We will need the following lemmas. 

Lemma 2.1. Let cp e C^(R). 
(i) //m = 0,1,2,..., then 



a=—2m 



tr' ,...^.. „ ^-f2™W.^ „ (-irm! 



(2.1) a.c. I '".' . ^ <^{t)dt = Cm,i<^^'^"'\0), c 



mA 



T[a/2) (2m)! 



(ii) If m = 1,2, .. ., then 

(2-2) «-c- hFr7^vit)dt = c^^^n——±idt 

a=l~2m J i (a/2) J t 

R R 

(2.3) = -Cm,2 [v^''^\t)l0g\t\dt 



where Cm,2 = (r(l/2— 'm)(2m — 1)!) ^ and the integral on the right-hand 
side of (2.2) is understood in the principal value sense. 

Proof. Both statements summarize known facts from [27, Chapter 1, 
Sec. 3]. For instance, (ii) can be proved as follows. Using the equality 

\tr' = ^. ^^^^ ,, i\tr'-' sgntr--'\ 

T{a + 2m — 1) 
we write the left-hand side of (2.2) in the form 

- a.c. ^, ^J^""],,.. ,,, (ir^^'"'^sgnt,y(^-^)(t)). 
a=l-2m L [a + 2m —1)1 (tt/2) 

The latter yields the principal value integral 

1 f^i—^dt 

r(l/2-m)(2m- 1)! 7 t 

R 

which coincides with (2.3). D 

Lemma 2.2. Let n > 2, \h\ < 1. 

(i) The integral 

1 

(2.4) g^{h) = -^ f\t-hr-'{l-t'Y^~'y'dt, Rea>0, 

T{a/2) J 
-1 

extends as an entire function of a and this extension represents a C°° 
function of h uniformly in a E K for any compact subset K of the 
complex plane. 
(ii) Moreover, 

(2.5) a.c. gM = n{n-l)/2). 

a=6—n 
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Figure 1. The Euclidean case. 

The proof of this lemma is given in Appendix. 

3. The Euclidean case. Derivation of the Inversion 

Formula 

We recall that our aim is to reconstruct a C°° function / supported 
in the ball B = {x EM."' : \x\ < R} provided that the spherical means 

(M/)(e,t) = ^ / fi^-ta)da, (C,t) edBx M+, 

O-n-l J 

are known for all spheres centered on the boundary dB of B (Fig. 1). 
We introduce the "back-projection" operator P that sends a function 
-F(^,t) on dB X M_|_ to a function {PF){x) on B by the formula 

(3.1) (PF)ix) = j^J F{^,\x-^\)d^, xeB, 

dB 

where d^ stands for the surface element of dB and \dB\ denotes the 
area of dB. 

3.1. The case n > 2. Consider the following analytic family of oper- 
ators 

(3.2) (iv"/)(e,t) = y J — r^,,/2) f^y^ ^^' 

B 

{C,t) edB XR+, Rea>0. 
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Lemma 3.1. Let f be an infinitely differentiable function supported in 

B = {x eiSJ' : \x\ < R}. Then 

(3.3) a.c. (PiV"/)(x) = A„ / , ^^^l , dy, 

a=3-n J \X — t/|" ^ 

B 

(3.4) A„ = (2i?)2-"7r-i/2r(^/2). 

Proof. For i?e a > 0, changing the order of integration, we obtain 

{PNy){x) = J fiy)Kix,y)dy, 

B 

where 

(3.5) = p^^i^ / ^^^rlP'o- 

o-n-i y r(a/2) 

I 12 I 12 

X — y Lc — \y\ 

a = r, h = ——. r. 

\x — y\ 2R\x — y\ 

By the rotation invariance of the inner product, the integral in (3.5) is 
independent of a and can be written as 
1 

0"„-2 



r(a/2) 



|t_/,|"-l (1 _ e)(n-3)/2 ^^ ^ a^^29a{h); 



cf. (2.4). Note that \h\ < 1 — 6 for some 6 > because x and y belong 
to the support of / and the latter is separated from the boundary of 
B. Hence, Lemma 2.2 yields 

a.c. (PiV"/)(x)= ^^^^''""-^ / , ^^^l^ a.c. gM dy 

a=3-n <Jn-l J \X — 1/|" ^ a=3-n 

B 

= ^n I r^H^2 dy, A„ = (2i?)2-"7r-i/2 r(n/2) 
J \x — y\^ ^ 

B 

(to justify interchange of integration and analytic continuation, the 
reader may consult, e.g., [52, Lemma L17 ]). D 
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Let US obtain another representation of the analytic continuation of 
PN°'f, now, in terms of the spherical means of /. 

Lemma 3.2. Let f be an infinitely diff'erentiable function supported in 

B = {x eMJ" : \x\ < R}, 

.^.. r, I d (_i)K2-i]r((n-l)/2) 

(i) If n = 3,5, ... , then 



d^. 

t=\x-^\ 



(3.7) jx.cjPNy){x) = ^ j D-^r~\Mf){i,t)] 

dB 

(ii) If n = 4,6, ... , then 
(3.8) 

2R 

a.c. {PN^f){x) = ^ fd^ ftD^''r~'{Mf){U)] \og\t'-\x-^\'\dt. 

a=3-n TT K J J 

dB 

Proof. Passing to polar coordinates, we have 

2R 

(iV"/)(e, t) = ^„-i I ^y=^^ (M/)(e, r) r^-^r 





r(a/2) 





V,{r)dr, V^,(r) = ^W^-i(M/)(e,rV^). 



Since the support of / is separated from the boundary of B, there is 
an e > such that (P^{t) = when r ^ {e,4:R'^ — e). Hence, ^^{t) can 
be regarded as a function in C^(m) and we can write 



R 

Now, Lemma 2.1 yields the following equalities. 

For n = 3, 5, . . . : 

.. (N'^n(ff^ A n(-'^(f^^ a (-i)^"-^^/M(^-3)/2)! 

a=3-n '5 {n — 6)\ 

For n = 4, 6, . . . : 

a.c. {N^f){^,t) = 6n,2 Lf-^\r)\og\r-e\dr, 
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Jn,2 



r((3-n)/2)(n-3)!' 
Combining these formulas with the backprojection P and noting that 
operations a.c. and P commute, we obtain 



a.c. {PNV){x) = ^ I Vt'\\^ ' ^\') dt 



For n = 3, 5, . . . : 

a.c. ,^ .. , ,,^ , , „ „ 

a=3-n \dB 

' dB 

For n = 4, 6, . . . : 

4i?2 



ac. (PiV"/)(x) = ^ f di f vf-'\T) \og\T-\x-i\'\dT. 

:=3-n |C'-d| J J ^ 



a=3~n 

dB 



These formulas give the desired result. D 

Comparing different forms of the analytic continuation in Lemmas 
3.1 and 3.2, we obtain the following statement. 

Lemma 3.3. Let f be an infinitely diff'erentiable function supported in 
the unit ball B = {x e M.'' : \x\ < R}, D = ^^^ . Then 

(3.9) A„ / , ^'^^l , dy 

J \x — y\"- ^ 



' l^,\D-\t'^-\Mf){i,t)\ 



t=\x-i\ 



dB 

if n = 3,5, ..., 



2R 
5^ 



- j di j tD'^-^r-\Mf){i,t)]\og\e -\x-i\'\dt, 



ttR'" 

dB 

^ if n = 4,6,..., 
where A„ and 5„ are defined by (3.4) and (3.6), respectively. 

The left-hand side of (3.9) is a constant multiple of the Riesz poten- 
tial of order 2 defined by 

2,,,^, r(n/2-l) [ f{y)dy 



(3.10) {I'f){x) - , „/2 / I |„_2 

and satisfying 

fc=l ^^k 
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Thus, we arrive at the following result. 

Theorem 3.4. Let f be an infinitely diff'erentiable function supported 
in the unit ball B = {x E R"^ : \x\ < R}, D = ^^ . 
(i) If n = 3,5, ... , then 

(3.12) fix) = <i A f D--^r-\Mf){i,t)] di, 

J t=\x-i\ 

dB 

"^"•^ " ART{n/2) • 
(ii) If n = 4,6, ... , then 

2R 

(3.13) f{x)=d^,2^ j di j tD--^r-\Mf){i,t)] \og\e-\x-i\^\dt, 

dB 

(_l)«/2-l7]--"/2 



d. 



n,2 



2R{n/2-l)\ ■ 

3.2. The case n = 2. Let P be the open disk in R^ of radius i? centered 
at the origin. In this section, for the sake of completeness, we reproduce 
(with minor changes) the argument from [18], keeping in mind that the 
Riesz potential of order 2 in the previous section is substituted by the 
logarithmic potential 

(3.14) ^j)^^) = l.J f^y)iog\x-y\dy, 

V 

satisfying AIJ = f. 

The following statement is a substitute for Lemma 2.2. 

Lemma 3.5. Let -I < h < I, a E S^. Then 



(3.15) ^, = / log|^-a-/i|rf^ = -27rlog2. 

SI 

Proof. Owing to rotational invariance, we can write 

(3.16) ,. = 2|^-^=_*. 

This integral is known; see, e.g., [17, p. 296], [18].^ D 



A more general integral was evaluated in [7, Lemma 6.1]. 
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Lemma 3.6. Let f be a C°° function supported in V. Then 

2R 

(3.17) (IJ)(x) = ^JJ{Mf){^,t)log\t'-\x-^f\tdtd^+Cf, 

&D 

c/ = — ^ / f{y) dy. 

V 

Proof. Let 

{NJ){U) = J f {y)\og\t'-\y-^f\ dy. 

V 

Changing the order of integration and making use of (3.15) with 

I |2 I |2 

X — y \x\ — \y\ 

\x — y\ 2n\x — y\ 

we obtain 

{PNJ){x)= J fiy)K{x,y)dy, 

V 

where 

k*{x,y) = ^;^ hog\\x - ^f - \y - ^f\ d^ 

dV 

— / log I \x\ — Ivl — 2f ■ (x — y)\ df 

dV 

= — {log{2R \x-y\) + log \h-0-a\)de 
2% J 

n 

= log(2_R |a; — yl) + — ^ = logi? + log |a; — i/|. 

ZTT 

This gives 

(3.18) {PNJ){x) = J f{y) log \x -y\dy + log R J f{y) dy. 

V V 

On the other hand, {PN^:f){x) can be expressed in terms of the spher- 
ical means. Indeed, passing to polar coordinates, we have 

2R 

(NJ) (e, t) = 2nj{Mf) (e,r) log \r' - t^\ rdr 
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and therefore, 

2R 

(3.19) {PNJ){x) = ^ j j {Mf) {^,r)\ogy-\x-^f\ rdrdC 

dV 

Comparing (3.18) and (3.19), we arrive at (3.17). D 

Lemma 3.6 allows us complete Theorem 3.4 in the following way. 

Theorem 3.7. Let f be an infinitely diff'erentiable function supported 
in the disk V = {x eM."^ : \x\ < R}. Then 

2R 



(3.20) f{x) = A\^-^l l{Mf){^,t)\og\t'-\x-^\'\tdtd^. 



dV 



Formula (3.20) can be formally obtained from (3.13) by setting n = 2. 
It coincides with formula (1.4) in [18]. 

3.3. Modified inversion formulas. We can replace / by A/ in (3.9) 
Since / is smooth and supp/ is separated from the boundary of B, 
then PAf = — /. Furthermore, since u{x,t) = {Mf){x,t) satisfies the 
Darboux equation 

n — 1 

Du = Au — Utt Ut = 

and A commutes with rotations and translations, then 

{MAf){x,t) = {AMf){x,t) = L[{Mf){x,-)m, Vx ER^, t> 0, 

where 

d"^ n — 1 d 

^ d^^ t di' 
see, e.g., [29, p. 17]. This reasoning and its obvious analogue for n = 2 
give the following modifications of inversion formulas (3.12), (3.13), 
and (3.20) with the same constant factors: 
(i) If n = 3, 5, . . . , then 






(3.21) fix) = <i j D''-'r-\LMf){i,t)] 

dB 

(ii) If n = 2,4,6,..., then 

2R 

(3.22) f{x) = dn,2 j dijtD--^[t--\LMf){i,t)] \og\e-\x-i\^\dt. 



dB 
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Formula (3.21) agrees with [54, formula (3.8)]. 

4. Spherical means and EPD equations 

Consider the Cauchy problem for the Euler-Poisson-Darboux equa- 
tion: 

n + 2a — 1 

(4.1) DaU = Au-Utt Ut = 0, 

V 

(4.2) u{x,0) = f{x), Mi(x,0) = 0. 

As in the previous section, we assume that / is a smooth function with 
compact support in the ball B = {x EM."' : \x\ < R}. If a > (1 — n)/2, 
then (4.1)-(4.2) has a unique solution 

(4.3) M(x,t) = (M"/)(x,t), xeR'', tGM+, 
where M°'f is defined as analytic continuation of the integral 

(MV)(a:,t) = ^-^t!|^ l(l-\y\Y-'f{x-ty)dy, Rea>0; 

ls/l<i 
see [8] for details. If a = 0, then M"/ = a.c.M°'f represents the 
spherical mean (1.1). 

Consider the following problem: 

Given the trace u (^, t) of the solution of (4-1) - (4-2) for all (^, t) G 
OB X M+, reconstruct f{x). 

To solve this problem we need some facts from fractional calculus; 
see, e.g., [55, Sec. 18.1] or [17, Sec. 9.6]. For i?ea > and r^ > -1/2, 
the Erdelyi-Kober fractional integral of a function (/? on R+ is defined 
by 

(4.4) (/»(t) = -p-— / it' - r^r'^r'^^Mr) dr, t > 0. 

i (a) Jo 

In our case it suffices to assume that (p is infinitely smooth and sup- 
ported away from the origin. Then I°:ip extends as an entire function 

of a and r], so that I^(p = cp, (/^) p = I^+a^-, 

Assuming x = ^ E dB and passing to polar coordinates, we obtain 

(4.5) u^{t) = (M"/)(e,t) = ^ r (n/^ ^ ^^^^^^ ^^^ ' 
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where u^ (t) = u{^,t), ip^ (t) = {M f){^,t), r] = n/2 — 1. This gives 
(A «^ r(n/2) . . -1 r(n/2) 

Now, since v^^(t) = {Mf){C,,t) is known, we can use Theorem 3.4 to 
reconstruct / by the following formulas. 

Theorem 4.1. Let f be an infinitely diff'erentiable function supported 
in the unit ball B = {x E R"^ : \x\ < R}, D 
(i) If n = 3,5, ... , then 



2t dt 



fix) = d^,, A I D-^[t--\l^^^u^) (t)] 



t=\x-(\ 



d. 



dB 

(_l)("-l)/2^1-n/2 

"'^ " 4:RT{a + n/2) 



in) If n = 2, 4,6,..., then 



2R 



fix) = dn,2 ^Jd^Jt ^"~'r~'(/,-+>5) it)] log \t' - k - en dt. 



dB 

(_l)«/2-l 7r-n/2 



(i. 



n,2 



2RTia + n/2) ' 



The case a = (1 — n)/2 in this theorem gives explicit solution to the 
TAT problem (see Introduction) with constant speed c(a;) = 1. More- 
over, after / has been found, we can reconstruct uix,t) in the whole 
space by setting ■u(x,t) = (M"/)(x, t). The latter gives an explicit so- 
lution to the Cauchy problem for the generalized EPD equation (4.1) 
with initial data on the cylinder dB x M+. 

5. Spherical Means on S'^ 

The suggested method of analytic continuation enables us to study 
the spherical mean Radon transform Mf on arbitrary constant cur- 
vature space X. In this setting, supp/ C B, where i? is a geodesic 
ball centered at the origin, and the spherical means of / are evaluated 
over geodesic spheres, the centers of which are located on the bound- 
ary of B. In this section we consider the case, when X = 5*" is the 
n-dimensional sphere in IR"^^. 

Given a; G 5" and t G (-1, 1), let 

(5.1) iMf)ix,t)=^- '- / fiy)daiy) 

x-y=t 
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Figure 2. The spherical case. 



be the mean value of a function / G C°°(S'") over the planar section 
{yeS^:x-y = t] (Fig.2). 

Our aim is to reconstruct / under the following assumptions: 

(a) The support of / lies on the spherical cap 

(5.2) Be = {x e S"^ : X ■ e„+i > cos 6*}, 

(the geodesic ball of radius 6), where Cn+i = (0, . . . , 0, 1) is the north 
pole of ^" and 6 e (0, 7r/2] is fixed. 

(b) The mean values (5.1) are known for all a; = ^ G OBq and all 
t G (— 1, 1), where dBe is the boundary of Bq. 

This problem can be solved using the method of the previous section. 
Let us fix our notation. In the following S^ = {x & S^ : x„+i > 0} 
is the upper hemisphere, Bq denotes the unit ball in the hyperplane 
Xn+i = 0; S*"^^ stands for the boundary of i?o, which is also the bound- 



ary of S*?. For x G S*? we write 



X 



{x\^/l^ 



X 



l\2\ 



X 



(a;i,...,a;„,0) G Bq, 



so that 
f\x) dx 



fix', v^ 



\X' 



'\2\ 



1 + 



dx 



n+l 



dxi 



+ ...+ 



dx 



n+l 



dXr 



dx' 



5? 



Bo 



fix', Vl - k'P) 



dx' 



Bo 



\X' 



'|2 
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We introduce the backprojection operator P, that sends functions 
on dBg X (—1, 1) to functions on Bq by the formula 

(5.3) {PF){x) = -j^ [ F{^,^-x)d^, xeBe. 



e 



Here d^ and \dBg\ denote the surface element and the area of dBg, 
respectively. We denote by Bg the orthogonal projection of Bg onto 
the hyperplane Xn+i = 0. 

5.1. The case n > 2. Assuming (^,t) G dBg x (—1, 1) and Rea > 0, 
consider the following analytic family of operators 

(5.4) (iv-/)(e, t) = J '^ "/("/Jp fiy) dy- 

Be 

Lemma 5.1. Let f e C°^(5"), supp/ C Bg. Then 

r(n/2)(sin^)2-" f f{y')dy' 



1/2 / \^l _ „/|n-2' 



(5.5) a.c. (PN''f)(x 

a=3-n ' '^ ' Tl^l^ J \x' -y 

Be 

(5.6) fiy') = il-\yf)-'/'fiy',il-\yfy/'). 

Proof. For Rea > 0, changing the order of integration, we obtain^ 

{PN^f){x) = Jf{y)Ux,y), dy, K{x,y) = ^ J ^l^-|I^rfe 

Be dBe 

Since C, has the form C, = Cn+i cos 6 + u sin 6, u & S"^^, then 

\^- {x-y)\ = I (x„+i - yn+i) cos 9 + (x' - y') -u sin 6'| 

(5.7) = \h — uj ■ a\ \x' — y'\ sind, 

(5.8) /,= ^!±i^l!±icot^, a '^' " ^' 



X — y \ \x' — y'\ 



Hence, 



, , , (Ix'-y'l sin^)"-i [ \h-uj-a\^-^ , 
(5.9) K{x,y)=^-^ ^ '- / ' ^, , ' du. 



For the sake of convenience, we use some notations which mimic analogous 
expressions in the Euclidean case. 
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The integral in (5.9) is independent of i^ and can be written as 



1 

Cn-2 I 



r(«/2) 



t_/,|"-l (1 _ f)(n-3)/2 ^^ ^ ar,.29a{h); 



-1 



cf. (2.4). This gives 

/P-TH^ 7 f \ ^n-2(|a:'-i/'| sin^)"-i 
(5.10) ka{x,y) = 9a{h). 

O'n-l 

Let us show that \h\ < 1. We write 

X = e„+i cos 7 + M sin 7, y = e„+i cos 5 + f sin S, 

7,0 6(0,6'); u,v & b , /i ^ 



f' ~ y'l 



Then 

(cos 7 — cos 5)^ 



\h\' 



|-usin7 — f sin(5P 

(cos 7 — cos 5Y 



< 



sin^ 7 — 2(-u ■ v) sin 7 sin 5 + sin^ 5 

(cos 7 — cos 5Y (cos 7 — cos 5Y 

sin^ 7 — 2 sin 7 sin 5 + sin^ 6 (sin 7 — sin 5Y 

Without loss of generality, suppose that 7 < 5. Then 

, ~ , cos 7 — cos 5 7 + (5 „ 

\h\ < = tan < tant^, 

sin — sin 7 2 

and therefore, \h\ = \h\ cot 9 < 1. 
Since \h\ < 1, Lemma 2.2 yields 

fiy)dy f f{y')dy' 



a=3-n J \x' — y'y^ ^ J \x' — y I" ^ 

Be Be 

D 
As before, we need one more representation of a.c. {PN'^f){x), now 

a=3~n 

in terms of the spherical means {Mf){C,,t). 
Lemma 5.2. Let f e C°^(5"), supp/ C Bg. Then 



ac. {PNy){x) = J^ Ud/dtr-%Mm,t){i-er''~'] 

a=3-n (smyj" ^ J 

dBg 



t=^-x 
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if n = 3,5, ... , and 

a.c. {PN''f){x) = , /" , / d^ 

dBe 
1 

X I {d/dtr-'[{Mm,t){i-t'r/'-']\og\t-^-x\dt, 

cos 29 

if n = 4,6, . . . , where 5„ is defined by (3.6). 

Proof. For Re a > 0, by making use of the formula 

1 
(5.11) jf{y)a{i-y)dy = <y.,,^,ja{r){Mf){i,r){l-^f^-UT, 

S"- -1 

we have 

1 

{NVm, t) = ^^ jiMfm, r) \r - tr\l - r^)"/^-^ dr 

-1 



Mr + t)dr, v^5(r) = a„_i(M/)(e,r)(l-r2)' 



r(«/2) 

R 

Since / is smooth and its support is separated from the boundary dBg, 
then {Mf){^, r) is smooth in the r- variable uniformly in ^ and vanishes 
identically in the respective neighborhoods of r = ±1. Thus, we can 
invoke Lemma 2.1 which yields the following equalities. 
For n = 3, 5, . . . ; cos 2^ < t < 1: 



a=S—n ^ 



For n = 4,6, 



a.c. (iV"/)(e,t) = -^ / ^f~'\r) log\r-t\dr, 

a=d—n TT J 

cos 29 

6n being defined by (3.6). The above formulas mimic those in the proof 
of Lemma 3.2 and the result follows. D 

Lemmas 5.1 and 5.2 imply the following inversion result for the spher- 
ical means on S'^. In the statement below, A^i = df + . . . + d^ is the 
usual Laplace operator in the x'-variable. 
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Theorem 5.3. Let f E C~(5"), supp/ C Bq. Then 
(5.12) /(x) = ^^A.,/o(x',v/T^M^), d. ^ '^ 



sine "^^^ 'V I I 7. -n 2"-iW2-ir(n/2)' 
where foix) = fo{x', a/1 — |a;'P) /ias t/ie following form: 

foix) = - f {d/dtT-%Mf){i,t){i-er"-'] di 

dBe 

if n = 3,5, ... , and 
1 

f,{x) = ^ j di j {d/dtr~\Mf){i,t){i-eri^~']\og\t-i-x\dt 

dBe cos 26 

ifn = A,6,.... 

5.2. The case n = 2. We keep the notation of section 5.1. Let 

(5.13) ilj)ix) ^ {hf){x', Vl-kf) = ^j f{y) log \x' - y'\ dy. 

Be 

SO that 

(5.14) A.,(/J)(x) = (1 - |x'n-V2 /(:^) = /(:^)/a;3. 

Lemma 5.4. Lei f he a C^ function supported in Bg. Then 

1 

(5.15) (IJ)(x) = j^^ j JiMf)i^,r)\og\T-^-x\dTd^+Cf, 

dBe -1 

Cf = -^ [log '-^) I fiy)dy. 

Be 
Proof. Let 

{NJ){i,t) = j f{y)\og\i-y~t\dy, (C, t) G 95, x (-1, 1). 

Changing the order of integration, owing to (5.7), we obtain 
{PNJ){x)= I f{y)K{x,y)dy, 
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where 

k*(^^y) = Tjr^-i ^og\^-{x-y)\d^ 
\o^e\ J 

dBg 

= -— / [log \x' — y'\ + log sin6' + log \h — u ■ a\] duo 
2% J 

= log \x' — y'\ + log sin 9 -\ — -, 

271 

//»1 1 I I 7 1 

log \h — cu ■ a\ dcu = 2 / — -r^=^^ dt = — 27rlog2; 
J-l V 1 — t^ 

cf. Lemma 3.5. This gives 

(5.16) {PNJ){x) = 2n{IJ){x) + (log ^) Jf{y)dy. 

Be 

On the other hand, by (5.11), 

1 

(5.17) (PiVJ)(x) = -i- [ [{Mm,T)\og\T-^-x\dTdC. 

smO J J 

dBg -1 

Comparing (5.17) with (5.16), we obtain the result. D 

Lemma 5.4 allows us complete Theorem 5.3 in the following way. 

Theorem 5.5. Let f be an infinitely differentiate function supported 
in the spherical cap Be = {x E S'^ : x ■ e^ > cos 6}, 9 G (0, 7r/2]. Then 

1 

(5.18) /(x) = ^-^A., I y"(M/)(e,r) log |r-e-x|cir rfe 

a Be -1 

Formula (5.18) can be formally obtained from (5.12) by setting n = 2. 

6. The inverse problem for the EPD equation on S*" 
The Euler-Poisson-Darboux equation on S*" has the form 

(6.1) Dq-m = 5xU — -u^^i^ — (n — 1 + 2a) cotuu^ + a{n — 1 + a)u = 0. 

Here a; G 5*" is the space variable, u G (0, vr) is the time variable, 
Sx is the relevant Beltrami-Laplace operator. For the sake of simplic- 
ity, we restrict ourselves to the case Rea > —n/2. In this case the 
corresponding Cauchy problem 

(6.2) U^u = 0, u{x, 0) = f{x), u^{x, 0) = 0, 
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with / G C°°{S"') has a solution u{x,ijj) = (M"/)(x, cos cj), where 
{M°'f){x,t) is defined as analytic continuation of the integral 

(6.3) iM-f)i^,t) = -^^—^^^^^^ J {x.y-tr-'f{y)dy, 

x-y>t 

c„,« = 2""V""/2r(a + n/2)/r(a), Rea>0, t G (-1, 1); 

see [44], [53, p. 179], and references therein. 

Let Bg = {x E S"- : X ■ e„+i > cos 9} be the spherical cap of a fixed 
radius 6 G (0, 7r/2] and let dBg be the boundary of Bq. Our aim is to 
solve the following 

Inverse problem. Suppose that the values g{^,u) of the solution 
of (6.2) are known for all {^,uj) G BBq x (0,7r). Reconstruct the initial 
function f G C°°(S'"), provided that the support of f lies in Bq. 

This problem can be solved using the results of the previous section. 
Assuming Re a > 0, we pass to spherical polar coordinates and write 

(6.3) as 

1 

(M"/)(e, t) = ^^ "_^';,^lZln/2 /(- - tr-\Mf){^, r) (1 - ryl^-' dr. 

t 
Then we set 

F^{t) = iMf){u)ii-t'r^'-\ 

Ga,dt) = p. ^ ;,. (1 - t'r-'^-/'g{^, COS-^ t), 

1 (a + n/2) an-i 
and invoke Riemann-Liouville fractional integrals [55] 

1 

(6.4) {I-u){t) = -^ [(t - tr-^u(T) dr, Rea>0. 

r(a) J 

t 
Thus, if -Re a > 0, then 

(6.5) (rF^)(t) = G.,^(t). 

Since / is infinitely differentiable and the support of / is separated 
from the boundary dBg, then F^ is infinitely differentiable on (—1, 1) 
uniformly in C, and suppF^ does not meet the endpoints ±1. It follows 
that (6.5) extends by analyticity to all complex a, and we have 

(M/)(e,t) = (l-t2)i-'^/2(/I"G«,^)(t) 

where /Z" is understood in the sense of analytic continuation. Now 
Theorem 5.3 yields the following explicit solution of our inverse prob- 
lem. 
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Theorem 6.1. Let Be be the spherical cap on 5" of radius 6 G (0, n/2], 

2l-a„n/2 

Ga,di) = rt ^ /o^ (1 - ^')""'^"^'^(e. COS"' ^)' ^e a > -n/2, 

where g is a given function on OBq x (0, vr). // the initial function f in 
the Cauchy problem (6.2) is infinitely differentiable and the support of 
f lies in the interior of Bq, then f can be reconstructed by the formula 

(6.6) f{x) = -^^A^.f,{x\^l^W\^), d^^ ^ '' 



sin^ -^^v >v 11/, n 2"-iW2-ir(n/2)^ 



where fo{x) = fo{x', a/1 — |a;'P) has the following form: 



d^ 



fo{x) = - j {d/dtr-\lZ^G^,^){t)] 

dBg 

if n = 3,5, ... , and 

1 

fo{x) = 1 J d^ J id/dtr-^[ir''G^,^){t)] log \t-^-x\dt 

dBg COS 26* 

z/n = 2,4,6,.... 

We recall that the case a = (1 — n)/2 in this theorem gives a solution 
of the relevant inverse problem for the wave equation on S*" in the 
framework of the formal spherical TAT model. 

7. Spherical Means in the Hyperbolic Space 

Most of the facts listed below can be found in [60]. Let E"'-*^, n > 2, 
be the real pseudo-Euclidean space of points x = {xi, . . . ,Xn+i) with 
the inner product 

(7.1) [x, y] = -xiyi x^yn + Xn+iyn+l- 

The hyperbolic space H" is interpreted as the "upper" sheet of the 
two-sheeted hyperboloid 

(7.2) H" = {x G E"'^ : [x, x] = 1, Xn+i > 0}. 

The hyperbolic coordinates of a point x = {xi, . . . ,Xn+i) G H" are 
defined by 

Xi = sinh r sin Wn-i ... sin ^2 sin Ui , 
X2 = sinh r sin Un-i . . . sin UJ2 cos Ui, 

(7.3) \ 

Xn = sinhr cosC(;„_i, 
Xn+i = coshr, 
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Figure 3. The hyperbolic case. 



where 

< Wi < 27r; <cjj < n, 1 < j < n - I; < r < oo. 

By (7.3), each x G H" can be represented as 

(7.4) X = u sinhr + e„+i coshr = {u sinhr, coshr) 

where a; is a point of the unit sphere S*""^ in M" with Euler angles 
Ui, . . . ,Un-i- We regard M" as the hyperplane x„+i = in E'^'^. The 
invariant (with respect to hyperbolic motions) measure dx in H" is 
given by dx = sinh"~^ r dcodr, where dco is the surface element of S^^^. 
The geodesic distance between points x and y in H" is defined by 

dist(x, y) = cosh~"'^[a;, y] (i.e., cosh dist{x,y) = [x,y]). 

Given a; G H" and t > 1, let 

(t2 - l)(l-")/2 



(7.5) 



iMf)ix,t) 



O-n-l 



f\y) d(r{y) 



[x,y]=t 



be the mean value of / over the planar section {y G H" : [x, y] = t}. 

As before, our aim is to reconstruct a function / G C°^(lHI") under 
the following assumptions: 
(a) The support of / lies in the geodesic ball (Fig. 3) 

B = {x eif : dist(a;,e„+i) < R} = {x eif : Xn+i < coshi?}, 

where e^+i = (0, . . . , 0, 1) is the origin of H" and i? > is fixed. 
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(b) The mean values (7.5) are known for all a; = ^ G dB and all t > 1, 
where dB is the boundary of B. 

Owing to (7.4), we can write a; G H" as 

x= (x, v^l+ |x'|2), x' = (a;i,...,a;n,0) G M", 
so that 



(7.6) f{x)dx= f{x',^/T+\^)p{x')dx', p{x')- 



l + 2|x'|2 

7i2~' 



l + \x 

We introduce the "back-projection" operator P that sends functions 
on dB X (1, oo) to functions on B by the formula 

(7.7) iPF){x) = T^J m. [e, x]) daiO, xeB. 

dB 

Let 

B = {x eR^ : \x'\ < sinhi?} 

be the orthogonal projection of B onto the hyperplane x„+i = 0. If 
^ = e„+i cosh_R + u sinhi?, u G 5*""^ and x = {x',Xn+i) G B, then 



l^jx] = a/ 1 + |x'|^ coshi? — (x' ■ u) sinhi?, 

and {PF){x) is actually a function of x' G B. We denote this function 
by (PF)(x'). 

7.1. The case n > 2. Let, as above, C, G (95, t > 1. Consider the 
analytic family of operators 

(7.8) ^Ny){U) = p-^^^^j^fiy)dy, Rea>0, 

B 

Lemma 7.1. If f e C^im"), supp/ C B, then 

(7.9) ... ,PA^"/)(.0 ^ ^'"/^"-f^'"" / ^M^ ..'. 

«=3— n TT ' J X — y \ 

B 



(7.10) /(^') = /(^', yiTb?) W Y 



2|f 



/|2 



b? 



Proof. For i?ea > 0, changing the order of integration, we obtain 
{PNy){x) = jf{y)K{x,y)dy, k^{x,y) = ^J ^Ml^lRlldaii). 

B dB 



SPHERICAL MEANS 25 

Since ^ has the form ^ = e„+i coshi? + u sinhi?, u G S""^^, then 

I [^, ix-y)]\ = I (a;„+i - y„+i) cosh R - (x' - y') ■ u sinh R\ 
= \h — u ■ a\ \x' — y'\ sinhi?, 

(7.11) /, = ^!^±i^^cothi?, a = ,""' ~ y\ . 

\x' — y'\ \x' — y'\ 

Hence, 

, , , (b'-y' I sinh /?)"-! /• I/i-wctI"-! , 

ka(x,y) = / ; — duj 

^ ' (Tn-i J r(a/2) 

an-2{\x' -y'\s\nhRY-^ 
= 9a{h), cf. (5.10). 

If \h\ < 1, we can apply Lemma 2.2 and write the integral over B as 
that over 5 C M"". This will give the result. 

It remains to show that \h\ < 1. By symmetry we may suppose that 
\y'\ < \x'\ , which we shall do from now on. Let 

a = \y'\, b = \x'\, bo = sinhi?. 

Since \x' — y'\ > \x'\ — \y'\, then 



h< fa [b) COthi?, fa (6) = 7 . 

b — a 
For a fixed, the function fa (b) is increasing in (a, oo), because 
D{a,b) 



fUb) = - ^^7=, D{a,b) = V{l + b^){l + a^)-l-ab>0. 

[b — a) vl + o 

Hence, h < /^(sinhi?) cothi?. The right-hand side of this inequality is 
less than 1. Indeed, setting bo = sinhi?, we have 

bo- a bo 

^ i + 6|^v2l±ZHl±M < 1, 

(bo -a) bo 



^l + bl-^{l + a^){l + bl) < {bo - a) bo, 
^0 < D{a,bo) . 
This completes the proof. D 
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Lemma 7.2. Let f G C°°(h"), supp/ C B. Then 



dB 



t=[i,x] 



if n = 3,5, ... , and 



dB 

cosh 2R 

X J id/dtr~'[{Mm,t){t'-ir/'-']iog\t-[^,x]\dt, 
1 

if n = 4,6, . . . , where 6n is defined by (3.6). 
Proof. For Re a > 0, by making use of the formula 

oo 

(7.12) j f{y) a{[^, y]) dy = a„_i J a{r){Mf){^, r) (r^ - If'^-' dr, 

M" 1 

we have 

oo 



^^iT + t)dr, v^^(r) = a„^i(M/)(e,r)(r2-l): 



r(a/2) 

R 

Since / is smooth and the support of / is separated from the boundary 
dB, then (M/)(^,r) is smooth in the r- variable uniformly in C, and 
vanishes identically in the respective neighborhood of r = 1. Thus, 
Lemma 2.1 yields the following equalities: 
For n = 3, 5, . . . : 

a=S—n ^ 



For n = 4, 6, . . . : 



cosh 2R 



a.c. (iV"/)(e,t) = -^ / v^f-2)(r) log |r-t| dr, 

1 

5n being the constant from (3.6). Now the result follows; cf. Lemmas 
3.2 and 5.2. D 
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Lemmas 7.1 and 7.2 imply the following inversion result for the spher- 
ical means on H". We recall that A^./ denotes the usual Laplace operator 
in the x'-variable. 

Theorem 7.3. Let n > 2. An infinitely diff'erentiable function f sup- 
ported in the geodesic ball B = {x E M"- : dist(x,e„+i) < R}, can be 
reconstructed from its spherical means {M f) (^, r), (,^, t) G (9-Bx (1, oo), 
by the formula 






|x|sinhi? -^^v 'V 11/' n 2"-iW2-ir(n/2)' 

where \x\ = a/I^'P + x'^_^_i and fo{x) = fo{x', a/1 + |x'P) has the fol- 
lowing form: 

fo{x) = - j{d/dtr-%Mf){i,t) {e - i)"/2-i] 

dB 

if n = 3,5, . . . , and 

cosh 2R 

f,{x) = ^jdi I {d/dtr~'[{Mf){u){t'-ir^'-']iog\t-[^,x]\dt 

dB 1 

ifn = A,Q,.... 

7.2. The case n = 2. The argument follows Section 5.2 almost verba- 
tim. Let 

(7.13) {IJ)ix) = i- y fiy) log \x' - y'\ dy, 

B 

SO that 

(7.14) A,,(/J)(x) = fix') = \x\ f{x)/xs. 
Lemma 7.4. /// be a C^ function supported in B, then 



(7.15) {hf){x): 



oo 
L. J J (Mf) {^,T)l0g\T - [^,x]\dTd^ + Cf, 



\dB 

dB 1 



Cf = -^[\og'-^) j f{y)dy. 



B 



Proof. Let 

(iVJ) (e,t) = j f{y)\og\[i,y]-t\dy, {i,t) e dB x (l,oo). 



B 
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Changing the order of integration, we obtain 

{PNJ){x)= [ f{y)K{x,y)dy, 



k^{x, y = log \x' — y'\ + log sinh R — log 2 (see the proof of Lemma 5.4) 
This gives 

(7.16) (PiVJ)(x) =27r(/J)(a:)+ (log ^^) lfiy)dy. 

B 

On the other hand, by (7.12), 

oo 

(7.17) iPNJ)ix) = — ^ j JiMm,T)\og\r-[^,x]\drd^. 

dB 1 

Comparing (7.17) with (7.16), we obtain (7.15). D 

Owing to (7.14), Lemma 7.4 allows us complete Theorem 7.3 as 
follows. 

Theorem 7.5. An infinitely diff'erentiable function f supported in the 
geodesic ball B = {x E B^ : dist(a;, 63) < R} can be reconstructed from 
its spherical means (Mf) {^,t), {^,t) G dB x (l,oo), by the formula 



(7.18) / (x) = ^-^p|-^ A., yy(M/)(e,r) log |r-[e,x] I rfrrfe 

dB 1 

Remark 7.6. As in Section 6, Theorems 7.3 and 7.5 can be applied to 
solution of inverse problems for the EPD equation in the hyperbolic 
space. The reasoning follows the same lines as before. We leave it to 
the interested reader. 



8. Appendix: Proof of Lemma 2.2 

It is convenient to split the proof in two parts. 
(i) We recall the notation 

1 

(8.1) ^„(/i) = — I— [\t-h\^-\l-ef^~^^/^dt, Rea>0, 
1 [a/ 2} J 
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where n > 2 and \h\ < 1 — 6, 6 > 0. Changing variables t = 2r — 1, 
h = 2C, — 1, we write ga{h) = ^^((l + h)/2), where 

1 

na+n—3 f 


(8.2) = f/,(0 + f/,(l-0, S/2<^<l-S/2, 



The last integral expresses through the Gauss hypergeometric function 
so that UaiO = a^{a) b{a) Ca{0^ where 

a^ia) = 2"+-3^(n-3)/2+.r((n - l)/2), 6(«) = ^^, 

1 „ /'^ — 13 — nn — 1 \ 

^"^^^ " r(« + (n-l)/2) V^~' ^~' ~^~ + "' ^J ' 
see, e.g., [51, 2.2.6(1)]. Owing to [16, 2.1.6], Cq(0 extends as an entire 
function of a, which is represented by an absolutely convergent power 
series. Since S/2 < ^ < 1—6/2, this series converges uniformly ina & K 
for any compact subset K of the complex plane. Furthermore, 0^(0;) 
is also an entire function and b{a) is meromorphic with the only poles 
— 1, —3, —5, . . .. Since g^ is an even function, i.e., ga{h) = ga{—h), these 
poles are eventually removable. Hence, Ga{C,) extends to all complex a 
as an entire function of a and this extension represents a C°° function 
of C, uniformly in a; G -ft'. This gives the desired result for ga{h). 
(ii) To compute analytic continuation oi ga a.t a = 3 — n, first, we 
assume 1/2 < Re a < 1 and |J?72a;| < 1 and represent Ga{0 ^^ ^ 
Mellin convolution 



00 

la+n— 3 



na+n-i r t P\ dr 

(8-3) ^«(^)=f(^/(0> f\0 = jh{r)h[^^- 


where 

r"+("-3)/2(l-r)("-3)/2 ifO<r<l, 
h{'r)={ h{T) = \l-rr\ 

0, ifl<r<oo. 

The Mellin transforms fj{s) = L fj{T)r^^^dT {j = 1,2) are evaluated 

as 

~ r(g + a + (n-3)/2)r((n-l)/2) 3-n 
/i s = ^;^ ■ , Res>^- Rea, 

T{s + a + n — 2) 2 
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T[s + a) r(l — s) 

By applying the convolution theorem and the relevant Mellin inversion 
formula [58], we obtain 

K+ioo 

where /(s) = /i(s)/2(s)- The function /(s) has poles in the half- 
plane Res < K at the points s = —j and s = —j — a — {n — 3)/2, 
j = 0, 1,2, . . .. Since 1/2 < Rea < 1, all these poles are simple, and 
the Cauchy residue theorem yields 

;«) . r (!i^) r,a.) f tf r r(. .^(n 3)/2) 
V 2 y ^ -^^ Lr(a-j)r(a-j + n-2) 

r+^"'3)/2 / r(-a + (3-n)/2-j) r((n-l)/2 + j) 

+ r((n - l)/2 - j) V r((3 - n)/2 - j) + r(a + (n - l)/2 + j) 
Ultimately, we arrive at the following expression for Ga{0'- 

(8.4) GaiO = \iF(l-a,3-a-n;^^-a;A 

X ^fS — nn— In — I \ 

where 

r((n - l)/2) (-1)" r(3 - a - n) sin avr 



23-"-"r(a/2) r((5 - n)/2 - a) cos (a + n/2)7r' 



r 



((n-l)/2) ^"+("-3)/2 r(a) / cosn7r/2 



Ai 



2 23-"-"r(a/2) r(a + (n - l)/2) V" ^ cos (a + n/2)7r 

Case 1. Let n = 2m, m = 2,3, . . .. Then 

where 

, , F(l — a, 3 — a — 2m; 5/2 — a — m; ,^) 

^^^'"'* ~ (-l)™r(a + 2m - 2) 1(5/2 - a - m) ' 

cot(«7r/2) F(3/2 - m, m - 1/2, a + m - 1/2; 

'^^'"^ ~ ^3/2-a-m r(l _ ^l) !(« + m - 1/2) ' 

A simple computation yields a.c. G'a(.^) =r(m— 1/2) = r((n— 1)/2). 

«=3— 2m 
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Case 2. Let n = 2m + 1, m = 1, 2, . . .. Then 
(8-6) G^iO = I^ZlI^^ 'l^'^ [^^(^' ") + ^2(^; ")]' 

2^ a /m QQg iycxTX 1 1) 

where 

F(l — a, 2 — a — 2m; 2 — a — -m; ,^) 






(_l)m+i r(a - 1 + 2m) r(2 - a - m) 
^"+™-iF(l - m, m; a + m; 



r(l — a) rfa + m) 



Passiiiff to the hmit as a — )• 2 — 2m, we obtain a.c. GniOi = r(m) = 

a=2-2m ^^ ^ ^ 

r((n — 1)/2). This completes the proof. D 

Remark 8.1. The basic equahty (8.4) can be proved in a different way if 
we rearrange hypergeometric functions in (8.2) using known formulas. 
Specifically, the second term in (8.2) can be transformed by formulas 
(33), (6), and (21) from [16, Section 2.9]. This gives 

^(0 = 7i('^) -^l 1 — tt,3 — n — a; a; ^ 



B{0 = (e(l-0r^""'^/S(«)i^(«,a + n-2;^ + a;e), 



, . ^ T{a + (n - l)/2) T{a + (n - 3)/2) 
^'^"^ ~ r(a)r(a + n-2) 

T{a + (n - l)/2) r((3 - n)/2 - a) sin(n - l)/2)7r 



^'^"^ '' r((n-l)/2)r((3-n)/2) sin(n - l)/2 + a)7r' 

Owing to [16, 2.9(2)], 

5(0 = 7.(«) e^"-)/^^" F (^, ^^■^ + a-^y 

This gives (8.4). 

An alternative proof of (ii) The following alternative proof is 
instructive and leads to the same result. In fact, it suffices to prove the 
equality 

(8.7) ac. gM = ^i^ - l)/2) 

a=i—n 
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in the weak sense. Indeed, suppose that 



(8.8) a.c. iga,i^)= a.c. gaih)tlj{h)dh = T{{n~l)/2) Hj{h)dh 

a=3—n a=3—n J J 

R B 

for any C°° function ip with compact support in the interval (—1,1). 
Since, by Part (i), the analytic continuation of ga{h) represents a C°° 
function of h uniformly in a G -ft' for any compact subset K of the 
complex plane, then (use, e.g., [52, Lemma 1.17 ]) 

ac. {ga,ip) = ( ac. ga,-^) 

a=S—n a=S—n 

and (8.8) yields ( a.c. ga,"^) = r((n— l)/2) {l,ip). This implies (8.7). 

a=3—n 

Let US prove (8.8). We denote 

where (■)+ stands for zero when the expression in brackets is non- 
positive. We interpret these functions as ©'-distributions on R. Then 
ga is a convolution of pa with the compactly supported distribution u 
so that {ga,ip) = (pa(s), (w(t),^(s + 1))); see, [27, Ch. I, Sec. 4(2)]. 
If n is odd, n = 2m + 3, m = 0, 1, . . ., then (2.1) yields 

(7 \ 2™ 
— ) (cj(t),V^(s + t)) 
as J s=o 

= c„,i(a;,^(2™)) = c™,i([(l-t^)™](2™),^(t)) 

= (-l)"c„,i(2m)! = m! = r((n-l)/2). 

Let now n be even. Since the convolution is commutative, 

a.c. {ga,i^) = a.c. (a;(t), (pa(s),^(s + t))) = {u(t),a.c. {pa{s),ij{s + t))). 

If n = 2m + 2 (m = 1, 2, . . .), then, applying (2.2) and changing vari- 
ables, we have 

a.c. {ga,ip)= a.c. {ga,V) = Cm,2 [uj{t),P-v. as) 

a=3—n a=l—2rn ^ J S / 

R 
1 1 

= (-l)'"+^c„,2 U^^"'-^\h)q{h)dh={-irCm,2 [^Pih)q^^"'-'\h)dh, 



-1 
1 



r(l/2 -m) (2m- 1)! J {t-h)^/l-t^ 
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(interchange of the order of integration can be justified if we complete 
[—1, 1] to a closed contour and use [25, Section 7.1]). Now, q{h) is 
a polynomial with leading term vr/i^"^"^. This follows from the well 
known relation for Chebyshev polynomials 

<«'^) "-/iT^fRib^ '«"-'<"'• -'<'^-'- 

and the fact that the leading terms oiTn(h) and Un{h) are 2"'~^/i" and 
2"/i", respectively; see formulas 10.11(47), 10.11(22), and 10.11(23) in 
[16, vol. II)]. Hence integration by parts yields 



a.c. (^«,^) = (-l)'"c™2vr(2m-l)! / ^lj{h) dh = T{m+l/2) / ^{h) dh, 

a=3-n ' J J 

-1 -1 

where r(m + 1/2) = T{{n — l)/2). Thus, we are done. D 
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